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The checkerboard pattern in the differential conductance maps on underdoped cuprates appears 
when the scanning tunneling microscopy is placed above the O-sites in the outermost Cu02-plane. In 
this position the interference between tunneling paths through the apical ions above the neighboring 
Cu-sites leads to an asymmetric weighting of final states in the two antinodal regions of fe-space. 
The form of the asymmetry in the differential conductance spectra in the checkerboard pattern 
favors asymmetry in the localization length rather than a nematic displacement as the underlying 
origin. 

PACS numbers: 74.72.-h, 74.25. Jb, 74.62.Dh 



I. INTRODUCTION 

The large scale detailed maps of the tunneling den- 
sity of states (DOS) in underdoped cuprates obtained 
recently by scanning tunneling microscopy (STM) have 
aroused great interest This surface sensitive technique 
to date has been limited to strongly disordered BSCCO 
and Ca2- a; Na a ;Cu02Cl2 samples. An analysis of these 
DOS-maps found a local breaking of the square symmetry 
of the underlying lattice to form a checkerboard pattern. 
The presence of this checkerboard pattern in these two 
different underdoped cuprates has been interpreted as ev- 
idence for the presence of an intrinsic bond-centered elec- 
tronic glass with unidirectional domains as a key charac- 
teristic of the pseudogap phase in underdoped cuprates P 
A number of different interpretations o f this patte rn have 
also been suggested in the literature! 3 * 4 * 5 * 6 * 7 *^ In our 
study we reexamine possible origins for this pattern in 
light of very recent STM experiments.^ 

The new STM data taken in the superconducting state 
at low temperatures, show that the checkerboard pattern 
appears most prominently at voltages where the tunnel- 
ing processes are predominantly into antinodal regions 
in fc-space. At lower voltages the spatial pattern of the 
DOS-maps is quite different. The low energy pattern 
in the superconducting state was successfully analyzed 
as arising from interference generated by the scattering 
of propagating Bogoliubov quasiparticles in the presence 
of weak disorder.^ Kohsaka et al. found that a rapid 
change in the spatial pattern in the DOS-maps occurred 
at the tunneling voltage corresponding to the energy of 
nodal centered arcs of Bogoliubov quasiparticles at their 
endpoints on the diamond formed by the lines in fc-space 
connecting the antinodal points.^Thus the checkerboard 
pattern which appears at higher voltages is formed by the 
tunneling of electrons and holes into pseudogap state lo- 
cated in the antinodal regions of fc-space. The key ques- 
tion to be answered is the origin and interpretation of 
this checkerboard pattern. 



The checkerboard pattern is characterized by a local 
symmetry breaking which reduces the square C4 symme- 
try of the Cu-O-Cu bonds in a (CuO)4 square plaquette 
to C2 symmetry leading to a unidirectional pattern of 
domains with a glassy short range order on the length 
scale of 4a (a: lattice parameter). This short range or- 
der with local symmetry breaking has led to proposals 
that a static spi n-charg e stripe glass coexists with the 
superconductivity! 5 -* 6 -*^* Another set of proposal d 3 * 4 * 7 * are 
based on the existence of static fluctuations in an or- 
der parameter which breaks both translational and C4 
symmetry leading to strong scattering of the Bogoli- 
ubov quasiparticles consistent with the observations. At 
this point we should comment on the difficulties associ- 
ated with an explanation based on an intrinsic symmetry 
breaking. While it is true that the STM experiments were 
performed on the highly disordered cuprates, BSCCO 
and Ca2- a; Na a ;Cu02Cl2, because of their good surfaces. 
Not all underdoped cuprates are disordered, for exam- 
ple, the two members of the YBCO family, YBa2Cu40s 
and YBa2Cu 3 06.5-ortho-II, are well ordered and intrin- 
sically strongly underdoped with hole concentrations es- 
timated to be x — 0.14 and 0.1 respectively. Indeed 
samples of these cuprates are of sufficient quality to al- 
low the observation of quantum oscillations at high fields 
and low temperatures.^ Zero field experiments on these 
well ordered samples show no evidence of static broken 
translational symmetry in nuclear magnetic resonance 
(YBa 2 Cu40 8 )^ or neutron scattering (YBa 2 Cu30 6 . 5 - 
ortho-II ) experiments.^ This lead us to conclude that 
the translational symmetry breaking observed by STM 
in the BSCCO and Ca2-^Na :! ;Cu02Cl2 samples results 
from the strong disorder associated with random doping 
a short range spin liquid in these cuprates. 

In this paper we propose that the checkerboard pattern 
observed by STM is caused by the local disorder in the 
Cu02 planes. Our explanation goes back to the obser- 
vation by Chen, Rice and Zhang (CRZ)P^ that when the 
STM tip is centered above the planar O-sites, it couples 



to two tunneling paths through apical CI ions above the 
neighboring Cu-sites and as a consequence interference 
between these paths plays an important role. This in- 
terference leads to a different tunneling matrix elements 
for electrons into hole states which are bonding (anti- 
bonding) upon reflection about the central O-site of the 
Cu-O-Cu bondP^ CRZ considered the dilute limit and 
examined the localized states of a single hole bound to a 
Na + -acceptor in a Ca2- a: Na 2 ;Cu02Cl2 sample. The Na + - 
acceptor sits above the center of a (CuO)4 square pla- 
quette. They found a doubly degenerate acceptor bound 
state which in the presence of a local quadrupolc electric 
field splits into two states, each with only C2 bonding 
symmetry. However this explanation of the local symme- 
try reduction from C4 to C'2 symmetry is not compatible 
with the recent STM data, which shows that the checker- 
board pattern is formed by the tunneling of electrons and 
holes into pseudogap states located in the antinodal re- 
gions of fc-space. The lowest energy for a single hole is 
at the nodal regions of fc-space so that single hole ac- 
ceptor bound states are made up from combinations of 
hole states near to the nodal not antinodal regions of 
fc-space. Nonetheless their main conclusion that inter- 
ference of tunneling paths occurs when the STM tip is 
above an O-site, remains valid. 

In light of this new STM dataP^we concentrate our at- 
tention on the role of the antinodal regions of fc-space as 
origins of the checkerboard pattern. There are two such 
regions, one near (7r,0) and the other near (0,7r). The 
bonding-antibonding pattern on the Cu-O-Cu bonds dif- 
fers between the two antinodal regions. For quasiparticle 
near the (0,7r) antinodal region the Bloch phase factor is 
bonding on Cu-O-Cu in the cc-direction and anti-bonding 
in the y-direction, as shown in FiglTTb) . This pattern 
is reversed for the case of the (n,0) antinodal region. 
Thus when the STM tip lies above the planar O-sites 
in a (CuO)4 square plaquette the stronger signal comes 
from the antinodal region with a bonding phase factor 
and a local difference in the tunneling DOS between the 
two antinodal regions will cause a breaking of the C4 to 
C2 symmetry. This symmetry breaking occurs only when 
the STM tip lies above the O-sites in a (CuO)4 square 
plaquette but not when the tip is above Cu-sites. Typical 
checkerboard patterns display strong modulation above 
the planar O-sites but much weaker modulation is ob- 
served at the Cu sites. This strongly suggests that we 
should concentrate on possible mechanisms that locally 
differentiate the two antinodal regions since it is this dif- 
ference which will show up at the planar O-sites. 

In their recent paper, Kohsaka et al. discussed the 
conditions that enhance the checkerboard pattern in the 
tunneling maps. They find that the pattern emerges most 
clearly in Z-maps which plot the ratio of differential con- 
ductances g(r, V"), at opposite bias, i.e. 

Z{r,E^eV) = g(r,+V)/g(r,-V) (1) 

Further the checkerboard pattern is most pronounced 
when the voltage V is adjusted locally to the local pseu- 
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FIG. 1: (Color online) Top panel: schematic demonstration 
of the interference of the possible two tunneling paths from 
STM tip to the Cu-O-Cu hybrid state when the tip is at the 
midway between two CI atoms. Lower panel: the contribution 
of the antinodal kA,y = (0, tt) state to the two interference 
paths. The blue and red dots are Cu sites with the alternative 
sign along y-direction denoting the sign of the Bloch wave 
function of antinodal k,A, y quasiparticles as shown in Eq[§] 
(The signs for the d x 2_ y 2 symmetry of Cu-orbitals are not 
shown.) For the k^y/x quasiparticle contributions to O y / x of 
the two paths cancel each other, meanwhile the contributions 
to O x / y of the two paths add together. 



dogap value. This value varies strongly across the field 
of view in the STM measurements. This is typically 
40a x 40a and the typical length scale of the variation 
in the local value of the pseudogap is roughly 4a, similar 
to the length scale of coherent checkerboard patterns. 



In this paper we examine possible sources of local sym- 
metry breaking from C4 to C2 symmetry in the STM 
patterns in the presence of strong local disorder, combin- 
ing the CRZ theory for STM tunneling with the Yang, 
Rice and Zhang (YRZ) ^ phenomenological form for the 
single electron propagator in the underdoped pseudogap 
phase. 



3 



II. YRZ-PHENOMENOLOGICAL SINGLE 
PARTICLE PROPAGATOR 

YRZ based their propagator on a generalization of the 
form of the propagator for a lightly doped array of 2- 
leg Hubbard ladders derived by Konik et aZM They 
introduced a self-energy, T, R (k,u), x), which diverges at 
u> = on a surface spanned by elastic particle-particle 
umklapp scattering analogous to the behavior of the lad- 
der model. In the 2D square lattice this umklapp sur- 
face is a diamond connecting antinodal points on the 
Brillouin zone boundary. Note this umklapp surface ap- 
pears as the energy gap surface also in the case of wider 
Hubbard ladders with more than 2 legs^. YRZ took 
over the form of d-wave RVB gap function A R (k,x) = 
Ao(x)(cos k x — cos k v ) from the renormalized mean field 
theory of Zhang et ai™ Using their Gutzwiller renormal- 
ization factor gt(x) for longer range hopping, YRZ pro- 
posed a form for the coherent part of the single particle 
propagator G p (k,u>,x) in the pseudogap state, 



G p (fc,w,x) 



9t{x) 



e(fc)- 

i 

A|(fc,z) 
cj + e°(fc) 



fi - E R (k,u,x) 



(2) 



where e (fe) only includes the NN renormalized hopping 
contribution, meanwhile e(fc) includes both the NN and 
longer range renormalized hopping contributions. 

The YRZ propagator fits well the quasiparticle dis- 
persion E\. and their weight Z l k seen in angle resolved 
photo emiss ion (ARPES) experiments on the pseudogap 
phase) 16 ! 20 ! It has aslo been used successfully to inter- 
pret a number of other experiments on the cuprate super- 
conductor, i.e. Raman scattering^, optical propertied, 
and specific heatP^. Recently Yang et alW^ showed that it 
also accounts well for the particle-hole asymmetry in the 
quasiparticle properties opening up along the Fermi arcs 
reported by Yang et al. in the normal pseudogap phased 
Very recentlj^H the temperature and doping dependence 
of superconducting penetration depth was also shown to 
fit well to this scheme. In the antinodal regions of fc-space 
the energy gap is largest with a magnitude A R (kA,x), 
where kA is an antinodal wavevector: kA.x = (f, 0) and 
k-A,y = (0,7r). The value of Aq(x) is taken from the 
renormalized mean field theory of Zhang et al. for the 
RVB stated It is a strong function of the hole density, 
x, dropping linearly with increasing x from Ao at x = 
to at a critical hole concentration x — £ c (~ 0.2). 



III. RANDOM ELECTRIC FIELD IN A 
THOMAS-FERMI APPROXIMATION 

In the cuprates BSCCO and Ca 2 _ ;r Na. I ;Cu02Cl2 which 
have been studied in the STM experiments, the holes are 



introduced by randomly distributed acceptors situated 
above and below the CuC>2-planes. The random electric 
field 4>{r) associated with the acceptors generates a ran- 
dom component in the local hole density. The case of 
a slowly varying potential can be treated in a Thomas- 
Fermi approximation. The hole chemical potential /i/j is 
a constant in space. In the YRZ propagator the hole 
pocket is determined by the quasiparticle energy contour 
at the chemical potential Hh which lies above the mini- 
mum hole energy E\ located at the midway between the 
node and (7r/2,7r/2) points in hole notation as shown in 
the upper panel of Fig. [2] The hole density of states 
(DOS) rises almost linearly with increasing hole energy 
in the YRZ ansatz, in agreement with recent angle in- 
tegrated photoemission dataPS The spatial varying hole 
Fermi energy EySr) relative to the bottom of the hole 
band in a simple Thomas- Fermi ansatz is given by 



E h (r) = fih- <t>{r) 



(3) 



leading to a spatial variation in the hole density x(r), as 
shown in the lower panel of Fig. [2] This has two effects, 
first, a breathing of the area enclosed by the hole pockets 
and therefore the Fermi arcs observed by ARPES and 
by quasiparticle interference in STM. Second, it intro- 
duces a local variation in the RVB gap A R (ktA, x) which 
as remarked above, is a strong function of the hole den- 
sity. The strong spatial variation in the RVB gap can 
be measured by the energy of the maximum in the DOS 
relative to the constant chemical p otential . That has 
been reported in several STM studies P 10 * 27 ! A breathing 
of the Fermi arc determined by quasiparticle interference 
in STM spectra, was recently reported by Wise et alW^ 



IV. INTERFERENCE IN TUNNELING 
SPECTRA WITH THE STM TIP ABOVE O-SITES 

The detailed STM maps showing the checkerboard pat- 
terns were taken on two compounds, Dy-doped BSCCO 
and Ca2- a: Na 2 ;Cu02Cl2. In BSCCO samples the out- 
most surface layer is a BiO plane with the Bi-ions located 
directly above the Cu-sites in the Cu0 2 -plane. Assum- 
ing the STM tips couple predominantly to the empty 
6s/6p Bi-orbitals, the analysis presented below will ap- 
ply also to BSCCO. We shall restrict ourselves to the 
former which has a simpler structure. At this point it 
is useful to recap the CRZ analysis, which discussed the 
role interference between neighboring tunneling paths in 
the STM maps. Meanwhile, we focus on the antinodal 
quasiparticles. With the STM tip centered at position 
r, the differential conductivity at T=0 and voltage V, is 
following Tersoff and HamannpHl 



dV ^ 

a.m 



| <m | a* | * > I 2 S(w - E m + E ) (4) 



where a\. a is the electron creation operator at position r, 
| m) are final eigenstates with energy E m , and u) = eV , 
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halfway between the two Cu-ions at i and i + t, the in- 
tegrated current will be sensitive to the relative phase of 
the hole states centered at the NN Cu -ions. Note the or- 
thogonality of the CI 3p 2 -states to the d x 2 _ y 2 -symmetry 
of the hole states centered on the Cu-site underneath 
leads to a dominant hybridization with the hole states 
centered on the 4 NN Cu-sitesP3 Finally, we arrive at 
the result for the tunneling currents when the STM tip 
lies halfway between NN CI sites, i.e. above the O-sites 
in the x/y oriented Cu-O-Cu bonds I°*/y(r + t/2,uj) 
and I Cu (r,ui) when the STM tip lies above a Cu site 

= e(uj-E m + E Q ) 



(in 



£(-i) t *(4 + xv + 4 + t+x') l*o) 



I Cu = Q(uj-E m + E ) 



(7) 



FIG. 2: (Color online) Top Panel: the low band from the 
single particle propagator Eq[2] in pseudogap state with the 
Fermi surface (hole pocket) projected on the xy plane. The 
bottom of the hole fermi sea is located at the midway between 
the node and (n /2,n/2). Lower panel: schematic demonstra- 
tion of the spatial variation of the hole Fermi energy Eh(r) 
due to the random electric field d>(r). 



| Wo) is the ground state with energy Eq. The outermost 
surface layer is composed of Cl-ions which sit directly 
above the Cu-ions in the topmost Cu02-layer, so that 
the STM tip couples primarily to their 3p z -states (see 
Fig 1). When the STM tip is scanned from above the 
Cl-ion at i to a NN CI- ion site, i + r [r = (±1,0) or 
(0, ±1)], we have 



Pom,. 



(*,Cl|r)^ M>(7 + <i + T,Cl|r)^ M+T) , 
^<i + r',Cu|i,Cl)4 +T/iCT 

t' 



(5) 



where (— l) r m is from the d x 



symmetry of the Cu- 



orbital, p cl i a is the electron creation operator at site i 

CI, and ct a is the electron creation operator for the d- 
p hybridized orbital centered on the Cu-site at i. The 
integrated current up to a positive voltage V (electron 
injection) is then 

I(r,u) = £|(m|^(-l)^[(i,Cl|r) C i +T> 



+ (i + r,C\\r)c\ +Tl+T \ |W 
Q(uj-E m + Eo) 



(0) 



is the Heaviside step function. CRZ pointed out that, 
when the tip is positioned above the O-ion at i+r/2 lying 



CRZ considered the case of a single hole bound to a Na + - 
acceptor which they showed may well have a doubly de- 
generate ground state. This degeneracy in turn splits 
in the presence of a quadrupole electric field, into states 
with reflection symmetries with respect to the x and y 
axes P x = ±1, P y = ±1 i.e. bonding and anti-bonding 
states at the O-sites. They showed that this lower C2 
symmetry shows up in the tunneling current. 

As we remarked earlier, the new STM data on un- 
derdoped Ca2- a; Na a ;Cu02Cl2 samples show that the 
checkerboard pattern is associated into antinodal states 
rather than the single hole bound states considered by 
CRZ. Note, however, quasiparticles in the antinodal re- 
gions states are composed of Bloch states near k,A, x = 
(it, 0) and fe^.j, = (0, tt), and these have similar reflection 
properties with opposite parities for reflections about the 
x- and y- axes. Focusing on the contributions from antin- 
odal quasiparticles, we can set 
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.(-ir*+cL,.(-i) 



(8) 



with the sign alternating pattern from h,A. y shown in 
Figjljb) . So that when the tip is above the O x / y site 
at 1 + t/2 



2 1) 1 [ C i+T',a + C i+-r'+T,a\ 

t' 

(9) 



it is obvious that for O x / y 
A,x/y quasiparticles vanishes, 

A,y/x 



From the factor 1 + (— l) r ^/» 
the contribution from the k 

meanwhile that from k,A :V / x is enhanced due to the in- 
terference effect, as schematically shown in Fig [l] As a 
result any local perturbation which breaks the symmetry 
between the two antinodal regions of fc-space will lead 
to a lower C'2 symmetry in the tunneling current. This 
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FIG. 3: (Color online) Typical checkerboard pattern revealed 
by R-map in Na-CCOC observed in one recent STM experi- 
ments RefP (reproduced from original paper) . The blue solid 
circles (red solid stars) show the position of the oxygen atoms 
on Cu-O-Cu bond along x/y directions [labelled by the com- 
bination of vertical lines (1-4) and horizon (a-i)]. 
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FIG. 4: (Color online) Differential tunneling conductance 
spectra taken on oxygen sites along the vertical lines (1,2,3,4) 
as indicated in Fig. [3] reproduced from Kohsaka et aZp^l 
Each subpanel shows the comparison of dl/dV between the 
two neighboring O^s (blue and cyan) and one O y (red) on 
each (Cu-0-Cu)4 palette (with the position shown in Figpl. 
Note that the spectra is normalized to have the same intensity 
at negative voltages to facilitate the comparison and demon- 
strate the lowering of the symmetry from C4 to C2. 



leads us then to look at various sources that can give rise 
to such local symmetry breaking in this highly random 
system. To proceed further it is useful to examine the na- 
ture of the anisotropy observed in the tunneling spectra 
more closely. 

V. MEASURED ANISOTROPIC STM 
TUNNELING SPECTRA. 

Kohsaka et a/P^ found that the checkerboard pattern 
emerges most clearly in Z-maps which plot the ratio of 
differential conductances g(r, V), at opposite bias, when 
the voltage V is chosen to be at the local value of the 
antinodal energy gap. In an earlier publication by the 
same group, individual spectra g(r } V), where plotted 
for a series of tip positions placed above O-sites in a 
small area with a pronounced checkerboard modulation P 
This is illustrated in the checkerboard modulation which 
shows up clearly in the R-map in Fig [3] which is repro- 
duced from their Fig.4BPThe R-ratio measures the ratio 
of the integrated currents I(r,V), at positive and nega- 
tive voltages 

R(r,E = eV) = I(r,+V)/I(r,-V) (10) 

The map reproduced in Fig. [3] shows the R-pattern taken 
at an energy E — 150 meV, with clear modulations when 
the tip is above the O-sites with a local C2 symmetry 
but almost no modulations are observed over the Cu- 
sites (denoted by black crosses) . The value of the R-ratio 
over the O-sites ranges from 0.9 (shown as white) to 0.5 
(black). 

The individual differential conductances taken above a 
series of O- and Cu- sites were also shown in their Fig. 



5B. We have extracted the individual tunneling spectra 
g(r, V) above the O-sites in this typical small area with 
a pronounced checkerboard pattern in the R-map (shown 
in Fig. 5B of Ref [2] ) . All spectra on these underdoped 
samples show a subst antia l asymmetry between negative 
and positive voltage;j- 30 * 31 l The negative voltage spec- 
tra, which correspond to electron extraction or hole in- 
jection, show little structure and has the larger differen- 
tial conductances. Since a hole can easily exchange posi- 
tions with a neighboring occupied site, while an electron 
can only hop onto the unoccupied neighboring sites, the 
tunneling processes that accompany an injected hole are 
much stronger than those for an electron, as shown by 
Anderson and Ong^ and by Randeria et al. Within 
the renormalized mean field theory (RMFT) of Zhang et 
aZP^, the asymmetry is largely due to the incoherent tun- 
neling process in the hole injection. In the approach we 
use here based on the YRZ normal state propagator, only 
coherent parts of the tunneling processes are included, 
and the particle-hole asymmetry appears at finite volt- 
ages which involves tunneling into quasi-particle states 
away from the chemical potential. The positive voltage 
spectra show a clear peak around an energy of order 100 
meV. Such a peak appears in the density of states calcu- 
lated for the YRZ propagator by Yang et al™ The YRZ 
propagator is an ansatz to describe coherent quasipar- 
ticles moving in a RVB background and so should agree 
better for positive voltage due to the weaker contribution 
of inelastic processes which accompany electron injection. 
In the subsequent discussion we will focus on the positive 
voltage spectra and normalize the spectra to be the same 
in the negative voltage region. In Fig. [4] we show a series 
of the differential conductances with the tip above O-sites 
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FIG. 5: (Color online) The density ol states on Cu and two 
kinds oxygen sites [O x / y ] from YRZ propagator Eq ' 
(panels a and b) and with (panels c and d) the ir 
factor (cosfca; — cosfc^) 2 . 



interference 

Panels (b, d) are the results in the 
presence of an intrinsic Ci asymmetry in the NN hopping in- 
tegrals along x and y direction t x / y /t = 1±0.1. The splitting 
of the peaks at positive voltage shows up clearly. The other 
parameters are as following: Ao = 0.4, t' = —0.3, t" = 0.2, 
doping x = 0.10 (all energies are of unit t — 300meV). Very 
small superconducting gaps are added to each band in pseudo- 
gap state, which does not change the RVB gap behavior (the 
two pronounced peaks at high energy ~ lOOmeV) on O x / y P 



in the x/y oriented Cu-O-Cu bonds. The difference from 
the lower C2-symmetry is clearly visible. 



VI. POSSIBLE SOURCES OF THE 
ANISOTROPY IN THE TUNNELING SPECTRA. 

We now discuss a number of physical effects which 
could cause an anisotropy in the tunneling spectra. 



Intrinsic C2 asymmetry due to nematic or 
Pomeranchuk order 



An asymmetry between the two antinodal regions near 
kA,x — (t, 0) and kA, y — (0, w) suggests that an in- 
trinsic instability which breaks the C4-symmetry may 
be present in these samples. The possibility of such an 
instability has been discussed extensively in the litera- 
ture. One set of these proposals focuses on a quantu m 
electronic liquid crystal with nematic order j32 P3 | 34 | 35 | j± 
second on a Pomeranchuk instability which splits the 
two Van Hove singularities! 36 * 3 ^ Recent neutron scatter- 
ing measurements of the low frequency spin fluctuations 
in YBCO samples have been interpreted as evidence for 
such intrinsic instabilities.^ In the presence of intrinsic 
asymmetry the NN hopping matrix elements in the x- 
and y- directions t x and t y will be slightly different. As- 



suming that the intrinsic asymmetry in the self energy 
in the YRZ propagator is small, the change of NN hop- 
ping integrals t x / y leads to a relative displacement of the 
reference energy for the RVB splitting at kA, x and kA, y - 
This in turn leads to a shift in the position of the energy 
gaps Afj(kA, x), relative to the constant chemical poten- 
tial. This is apparent in Fig. [5]which shows the predicted 
anisotropy arising from a uniform hopping asymmetry. 

The density of states, N(E), observed by putting the 
STM tip above O x / y and Cu sites takes the following 
form from Eq.(|7]) 

k 1 

N°*'y(E) ~ ^Z^cos 2 ^(cosfcz - cos k y ) 2 5(E l k - E) 
N Cu (E) ~ ^ Z l k (cos k x - cos k y ) 2 S(El ^ E) (11) 

i . k 

where (cosk x — cosk y ) 2 is the interference factor due to 
the hybridization with the hole states centered on the 4 
NN Cu-sites as discussed earlier. 

In a clean well ordered sample with intrinsic C2 asym- 
metry order large domains should occur but in the 
highly disordered Ca2_ a; Na 2 ;Cu02Cl2 samples there is a 
strong spatially varying electric field. The linear cou- 
pling between the orientation of the domains and the 
local quadrupole component of a random electric field 
acts as a random orienting field on the order. Since the 
symmetry of the order is Ising like, the random orient- 
ing field will act similarly to a random field in an Ising 
model. It is well known that the ground state of this 
Ising model consists of a random array of domains, con- 
sistent with the orientational domains observed in the 
checkerboard patterns. A word of caution, however, is in 
order. The local anisotropy in the STM spectra does not 
seem to be consistent with that predicted for this type 
local asymmetry order. The spectra in Fig. [4] show a C2 
-modulation pattern arising from an overall reduction in 
the positive voltage spectra, rather than a shift in the 
maxima that would signal the presence of this intrinsic 
C2 asymmetry in Fig (5] 

The theoretical DOS for the STM tip located above 
Cu- and O-sites shown in Fig|5] display considerably 
sharper structure in both cases, with and without the 
weighting factor (cos k x — cosk y ) 2 , than is evident in 
the experimental conductances in Fig|4j The theoret- 
ical DOS contain two prominent features. One is the 
approximately linear drop from negative to positive volt- 
ages leading to the minimum at ~ 60mcV. Interestingly 
this linear feature has been reported in integrated pho- 
toemission studies on BSSCO samples by Hashimoto et 
al.^as discussed recently by Yang et aZpSl. In the STM 
spectra there are signs of a minimum at positive voltages 
in the spectral 2h, 2f in Figj4] but at a smaller voltage and 
faint signs in the other spectra. The second feature is the 
very sharp peak at ~ lOOmeV in the theoretical spectra 
and the much broader peaks at around the same energy 
in the STM spectra. For both discrepancies a possible 
remedy could be the much stronger disorder broadening 



7 



in the local probe in the STM spectra as compared to 
spatially averaged spectra observed in photoemission. 

B. Random electric fields 

The random distribution of Na + -acceptors in 
Ca2_ a: Na 2 ;Cu02Cl2 samples generates a random electric 
field. The local quadrupole component of this field 
acting on a (CuO)4 plaquette will couple differently to 
the states at k^x and kA, y - The former is anti-bonding 
on the Cu-O-Cu bonds along the a;-axis and bonding 
on the Cu-O-Cu bonds along the y-axis while the latter 
has the opposite bonding pattern. These patterns cause 
opposite quadrupole charge distributions in the two 
states leading to a potential splitting of the antinodal 
states in the presence of an external electric quadrupole 
field. This potential splitting will generate STM spectra 
similar to those we found for hopping asymmetry and so 
it also differs from the anisotropy displayed in Fig. [4] 

C. Variations in the quasiparticle weight and in 
the localization length 

The tunneling DOS into quasiparticle states is 
weighted by the quasiparticle weight in the single par- 
ticle propagator. In the YRZ ansatz this weight is taken 
from the RMFT calculations of Zhang et aZPand is sim- 
ply proportional to the hole density, x. However the fact 
that little change in the energy gaps A^fcy^, x), is evi- 
dent in the spectra shown in Fig. [4] argues against such 
an explanation. Also the spatial variation in the hole 
density is longer ranged than the distance between O- 
sites on a (CuO)4. Another possibility is that there is 
substantial anisotropy in the disorder scattering in fc- 
space which causes a local anisotropy in the localiza- 
tion lengths. An anisotropic energy shift and localiza- 
tion length could result, if the potential scattering from 
the disorder is mainly between opposite antinodal points 
rather than between ifc^x and i:kA,y- An anisotropic 
energy shift, for example, could remove weight from the 
energy range ( 200meV) covered in the STM experiments. 
The somewhat longer range of the disorder potential is 
consistent with this suggestion. The umklapp correlated 
scattering wave vector for backscattering between the 



antinodal regions near ±kA, x (or ±kA, y ) is much smaller 
than the wave vector connecting the two antinodal re- 
gions, ~ (7r,7r)- The Fourier transforms of the checker- 
board maps show peaks arising from short range order 
in the form of 4a wide unidirectional electronic domains. 
The analysis presented here does not directly address the 
presence of short range order on a length scale > a. Sev- 
eral authors have suggested that the wave vector of this 
short range order is approximately equal to the wave vec- 
tor connecting parallel pieces of the bare band structure 
Fermi surface near the saddle points^ and also the wave 
vector connecting the turning points of the quasiparticle 
contours in the normal state YRZ propagators!^! The 
length scale of the domains at ~ 4a is roughly the same 
as the disorder potential length scale but the interpreta- 
tion of the short range order in the checkerboard patterns 
remains to be an open issue. 

VII. CONCLUSION 

In this paper we have examined several possible ori- 
gins for the local checkerboard pattern in the tunneling 
DOS measured over O x / y -sites. The presence of inter- 
ference between two tunneling paths centered above the 
neighboring Cu-sitcs leads to an anisotropic weighting 
of the two antinodal regions in fc-space. Thus the lo- 
cal C*2 pattern of the checkerboard can reflect an asym- 
metry between the two antinodal fc-space regions. An 
asymmetric hopping (t x , t y ) could arise from an intrinsic 
instability e.g. towards nematic or Pomeranchuk order. 
This would split the energies of the pseudogaps at ±kA, x 
and ±ftA,y However, an examination of the measured 
anisotropy in the differential conductance points rather 
towards an asymmetry in the average magnitude of the 
conductance, which possibly could result from a differ- 
ence in the localization lengths at the two antinodal sad- 
dle points in the quasiparticle dispersion. 
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